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1 INTRODUCTION 

Pulsars are high-velocity compact stars with transverse 
speeds ranging from a few tens of km s~* to over 1,000 km 
s _1 . These speeds significantly exceed those of their pro- 
genitor stars, which are typically of the order of 10 km s _1 
(Lyne & Lorimer 1994). A number of scenarios have been 
put forward to explain the high velocities of pulsars. Classi- 
cally, an asymmetry in the supernova (SN) explosion — of 
the order of 1% — is assumed to provide a kick that im- 
parts the high velocities observed today (Shklovskii 1969). 
Alternatively, the "rocket" model of Tademaru & Harrison 
(1975), in which the pulsar's magnetic-dipole configuration 
is offset from the progentitor star's centre, can provide birth 
velocities of ~ 100 km s - . The much higher actual veloc- 
ities observed and the requirement of this model for birth 
spin periods of Po ~ 1 nis have made it less preferred. How- 
ever, an interesting prediction of the rocket model is that 
the pulsars at birth are accelerated in the direction of their 
spin axes, thus implying a correlation between the angular 
and linear momenta. 



ABSTRACT 

The reported alignment between the projected spin-axes and proper motion direc- 
tions of pulsars is revisited in the light of new data from Jodrell Bank and Effelsberg. 
The present investigation uses 54 pulsars, the largest to date sample of pulsars with 
proper- motion and absolute polarisation, to study this effect. Our study has found 
strong evidence for pulsar spin-velocity alignment, excluding that those two vectors 
are completely uncorrected, with > 99% confidence. Although we cannot exclude 
the possibility of orthogonal spin-velocity configurations, comparison of the data with 
simulations shows that the scenario of aligned vectors is more likely than that of the 
orthogonal case. Moreover, we have determined the spread of velocities that a spin- 
aligned and spin-orthogonal distribution of kicks must have to produce the observed 
distribution of spin-velocity angle offsets. If the observed distribution of spin-velocity 
offset angles is the result of spin-aligned kicks, then we find that the distribution of 
kick- velocity directions must be broad with a v ~ 30°; if the orthogonal-kick scenario 
is assumed, then the velocity distribution is much narrower with a v < 10°. Finally, in 
contrast to previous studies, we have performed robustness tests on our data, in order 
to determine whether our conclusions are the result of a statistical and/or systematic 
bias. The conclusion of a correlation between the spin and velocity vectors is inde- 
pendent of a bias introduced by subsets in the total sample. Moreover, we estimate 
that the observed alignment is robust to within 10% systematic uncertainties on the 
determination of the spin-axis direction from polarisation data. 
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More recently, Spruit & Phinney (1998) and Cowsik 
(1998) proposed a mechanism whereby a number of impulses 
imparted onto the pulsar during core collapse provide not 
only the expected high birth velocities but also the short 
birth spin periods (Po ~ 10 ms). In this model, if only a 
single impulse of duration r is delivered during the collapse, 
then for t < Po this scenario leads to a spin axis that is 
preferentially orthogonal to the velocity vector (S _L v); 
however, if the pulsar receives a number of such kicks, this 
correlation is destroyed. Interestingly, if the the duration of 
these impulses is of the order of the pulsar's spin period or 
higher, then the spin axis is expected to be aligned with the 
velocity vector S || v); a similar conclusion was reached with 
toy-model simulations by Wang, Lai & Han (2007). 

Finally, since half of the isolated neutron stars are ex- 
pected to have come from binary systems, another mech- 
anism that has been proposed to explain the high pulsar 
velocities is the break-up of such binary system during the 
second SN (Gott, Gunn & Ostriker 1970; Bailes 1989). How- 
ever, break-up alone cannot account for observed velocities 
in excess of 1,000 km s - , so a significant fraction of the 
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observed velocity must also originate from the core collapse 
phase. At the moment, it seems that neither the simple con- 
servation of angular momentum nor the binary break-up are 
enough individually to impart the observed magnitude of 
pulsar velocities. Nevertheless, there are alternative physi- 
cal processes that have been recently suggested as the trigger 
mechanism for the presence of asymmetries in the SN ejecta 
(Lai et al. 2001; Janka et al. 2005). Those mechanisms can 
produce fast-moving pulsars without requiring short birth 
periods nor very strong magnetic fields. 

In addition to the origin of the magnitude of the pulsar 
velocities, another long-standing issue is the relative direc- 
tions of the pulsar velocities to the pulsar spin. As mentioned 
above, certain SN kick mechanisms predict a correlation be- 
tween the spin and velocity vectors of pulsars (e.g. Tade- 
maru & Harisson 1975; Spruit & Phinney 1998). In the case 
of pulsars whose kinematics and orientation are the result 
of the disruption of short-lived, close-binary systems, it has 
been suggested that the dominant kick component is due to 
the impulse received during the break-up of the system: it is 
therefore predominantly orthogonal to the spin axis (Colpi 
& Wasserman 2003). A weaker kick component is associated 
with a number of internal kicks during the core-collapse pro- 
cess, which tend to align the kick velocity with the spin axis 
(Spruit & Phinney 1998). The above authors present the 
case of the double neutron-star system, PSR B1913+16, as 
an example of the superposition of those two components: 
Kramer (1998) used the detection of geodetic precession in 
that system to measure the misalignment angle between the 
pulsar's spin and its orbital momentum (spin tilt); he con- 
sidered that as direct evidence for an asymmetric kick dur- 
ing the second SN. Later, Wex et al. (2000) concluded that, 
based on the orbital dynamics at the time of the second SN 
explosion, which they constrained using the system's proper- 
ties observed at present, the second SN kick must have been 
directed almost orthogonally to the progenitor's spin. Sim- 
ilar conclusions were later drawn by Willems et al. (2004), 
using and updated value of the spin tilt of PSR B1913+16, 
which confirmed that it is unlikely that the second SN kick 
was closely aligned with the progenitor's spin. 

The observational consequence of the above work is that 
fast pulsars (~ 500 km s _1 ) originating from close bina- 
ries that were disrupted by an asymmetric SN kick — but 
whose pre-SN orbital-momentum and spin vectors were all 
aligned — should have space velocities nearly orthogonal 
to their spin axes. On the other hand, pulsars with rela- 
tively low velocities (~ 100 km s _1 ) — i.e. those born in iso- 
lated SN events that imparted a number of weaker kicks on 
the proto-neutron star — could predominantly show spin- 
velocity alignment. 

This distinction between a population of slower pulsars 
and that of faster ones fits the model of Arzoumanian et 
al. (2002), who used a bimodal distribution with Gaussian 
components centred at 90 and 500 km s _1 to describe the 
distribution of pulsar birth velocities, based on the available 
sample of 435 isolated Galactic pulsars. However, the recent 
pulsar-population synthesis by Kuranov, Popov & Postnov 
(2009; hereafter KPP), who simulated the spins and veloci- 
ties of pulsars born in both binary systems and from isolated 
stars, suggests that the observed distribution of velocities is 
dominated by the SN kick mechanism, with the effect of bi- 
nary progenitors being negligible. This is indeed consistent 



with the observed single-mode distribution of 3D velocities, 
as was derived by Hobbs et al. (2005), from 233 pulsar proper 
motions. One of the consequences of the kick model used in 
KPP was that the origin of tightly aligned spin-velocity con- 
figurations (with alignment angles of S±Zv < 10°, where 
S± is the spin orientation) was mainly found for pulsars 
born in isolated SNe, having kick velocities of > 50 km s _1 . 
In contrast, binary progenitors resulted in misaligned con- 
figurations, covering the entire range of alignment angles. 

Ng & Romani (2007) performed detailed simulations of 
the kick dynamics during the core-collapse event — indepen- 
dently of any pre-existing systemic velocity (proper/orbital 
motion) — and concluded that fast-moving pulsars are more 
likely to show spin-velocity correlation than slower-moving 
ones. The parameters of their model were determined by fit- 
ting a sample of pulsars that included pulsars with (a) avail- 
able proper motions but no projected spin-axis orientation, 
(b) available proper motions and projected spin-axis orien- 
tations (e.g. from radio-polarimetry measurements), and fi- 
nally (c) available proper motions and knowledge of the 3D 
orientation of their spin axis (from the measured inclina- 
tion of the associated pulsar- wind tori). Although a pre-kick 
spin was required in those simulations in order to produce 
alignment, a short initial spin period (P < 20 ms) was not 
a strict prerequisite: the simulation produced also aligned 
spin-velocity vectors from long-period proto-neutron stars 
(P > 100 ms). Another interesting statistic that was borne 
out from that study was the standard deviation of the dis- 
tribution of the kick-spin alignment angles, a v , that led 
to the observed distribution of kick-spin offsets: the best- 
determined values were in the range a v « 10° — 20° . Similar 
results were derived by the aforementioned population syn- 
thesis of KPP; although it should be noted that comparison 
with pulsar data by Rankin (2007) led to a narrower distri- 
bution, i.e. a v ps 5°— 10°. 

This paper is organised as follows. In Section [2] we 
present the relevant, observational background information 
that needs to be considered in our study. Furthermore, in 
the same section, we describe the radio-polarisation obser- 
vations and the selected data that were used in this paper. In 
Section [3] we perform a statistical analysis on the selected 
sample and make some concluding remarks about the ob- 
served distribution of spin-velocity angles. In Section 4, we 
attempt to fit simple models of the birth-kick distribution 
to the observed distribution of angles; moreover, we pro- 
vide estimates of the width of the angle distributions, for 
the cases of aligned and orthogonal kicks. At the end of this 
section, we perform a systematics check for the robustness 
of our conclusions against various degrees of deviation from 
the published polarisation-angle values; in addition, we try 
to quantify the systematic influence of the inclusion of the 
Johnston et al. (2005) data — those being the most convinc- 
ing evidence for spin-velocity alignment — in our data set. 
Finally, in Section 5 we summarise our results. 



2 POLARISATION OBSERVATIONS 
2.1 Geometry 

Ideally, the determination of a pulsar's magnetic-pole ori- 
entation, having a north-through-east position angle (PAo) 



© 0000 RAS, MNRAS 000, 000-000 



Pulsar Spin-Velocity Alignment 3 




Figure 1. (a) Cross-section of the pulsar beam, showing the sky-projected orientation (radial, solid lines) and direction (arrows) of the 
magnetic field lines, and the trajectory of the observer's line of sight (dash-dotted line), as the beam sweeps past the field of view. The 
plane of polarisation of the normal mode of emission at each point of the line of sight is the plane that is perpendicular to the page 
and contains the corresponding trace of the magnetic field line. The orthogonal mode of emission occurs perpendicularly to that plane, 
(b) The characteristic S-shape of the position angle (PA) rotation across the pulse. According to the rotating vector model (RVM), 
as the pulsar beam sweeps the observer, our line-of-sight successively intersects the different magnetic field lines, corresponding to the 
different planes of linear polarisation, at different angles, hence producing the PA-profile shape shown, (c) 3D geometry of the relative 
configuration between the spin, magnetic and velocity axes. The beam geometry can be described by two angles: a, the angle between 
the spin and magnetic axes, and 0, the angle between the magnetic axes and the line-of-sight at the phase of minimum approach (<j>o)- 
The observer can only determine the direction of the sky-projected pulsar spin and velocity vectors, Si and v±, which are respectively 
defined by the position angles PAo and PA V , measured north-through-east. At <po, the plane containing the observer and the spin axis 
also contains the magnetic axis, thus allowing us to determine the spin-axis orientation. Finally, the opening angle of the emission beam 
is defined as p. 



on the sky, comes from a rotating-vector-model (RVM) fit 
to the characteristic 'S' shape of the radio-polarisation po- 
sition angle (PA) profile. Fig. []Ji,b shows schematically the 
geometric principle behind the rotating vector model, which 
was developed by Radhakrishnan & Cooke (1969). Such a 
model predicts that the point of inflexion in the PA profiles 
would correspond to PAo and it would therefore be straight- 
forward to determine PAo for any pulsar. However, there are 
few cases (e.g. the Vela pulsar, PSR J2043+2740) where the 
RVM swing can be clearly detected and then reliably fitted. 
Everett & Weisberg (2001) attempted fits to 70 pulsars' po- 
larisation profiles and could only reliably fit 10 cases. The 
main issue is that due to the narrow, conal radio emission, 
we can only obtain significant polarisation from pulsars over 
only a small fraction of their periods, whilst the rest of the 
magnetic-field-line rotation remains untraceable. 

As is shown in Fig. the RVM model describes the 
change of PA as a function of phase, <j>, using these param- 
eters: the angle between the spin and magnetic axes, a; the 
angle between the observer's LOS and the magnetic axis, /3; 
and the PA and (j> a t the minimum approach of the mag- 
netic pole to the observer's LOS, PAo and <f>o, respectively: 
i.e. tan(PA — PAo) = sinasin((jf> — </>o)/[sin(a + /3)cosa — 
cos(a + /3) sin a cos(0 — <j>o)]. For those pulsars where the 
RVM model can satisfactorily describe the observed PA pro- 
file, the relative orientation between the pulsar spin, mag- 
netic axis and the observer's line of sight (LOS) can be de- 
rived. 

Despite the limitations of RVM fitting, if one is only 



interested in determining PAo, there are alternative argu- 
ments that can help: e.g. one can use other indicators, like 
the pulse profile's centre of symmetry or the change of hand- 
edness of the circular-polarisation profile (see Section [2.1.11 
below). The latter feature, when it occurs near the profile's 
peak, has been associated with emission from near the mag- 
netic pole and viewing angles almost along the local mag- 
netic field (Rankin 1986). In the present work, when RVM 
fitting alone was deemed unreliable with our data, we also 
considered arguments based on pulsar symmetry, which took 
into account the distribution and size of the profile's peaks, 
as was done in Johnston et al. (2005, 2007). Such geomet- 
rical investigation benefited from available multi-frequency 
profiles for the pulsars in our sample, in order to account for 
a possible range in emission heights (Kramer et al. 1997). 
Overall, in most cases the determined PAo was a product of 
all the above considerations. 

It is also worth acknowledging that any of the above 
approaches of determining PAo may potentially lead to er- 
roneous results due to aberration/retardation (A/R) effects, 
which, although small, can potentially cause the inflexion 
point of the RVM swing to appear delayed with respect to 
the pulse profile's mid-point (Blaskiewicz et al. 1991; Dyks 
et al. 2004). In addition, A/R effects can shift the positions 
of conal components to earlier observed phases, thus mis- 
leading arguments based on profile symmetry. Clearly, ob- 
servations of rapidly-rotating pulsars will be more prone to 
the above effects. Recently, Krzeszowski et al. (2009) — and 
earlier Mitra & Li (2004) — investigated the presence of A/R 
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effects in pulsar profiles from multi-frequency observations. 
Amongst their conclusions was that the peak emission does 
not coincide with the fiducial phase containing the magnetic 
and spin axes and that, overall, the measured phase lags due 
to A/R were of the order of 1 deg. 

Upon measuring PAo, which is always measured relative 
to the orthogonal frame defined by our polarisation feeds, we 
have effectively determined an orthogonal reference frame 
on the sky plane that defines the spin-axis orientation. This 
becomes clear if we assume that the plane of polarisation 
of the measured signal (defined by PA) coincides with the 
plane defined by the magnetic-field lines and the magnetic 
axis, as is shown in Fig. []J. At the phase of minimum ap- 
proach of the observer's LOS to the magnetic pole (i.e. at 
PAo), that plane should also contain the spin axis. So, in 
that way, we can say that we have determined the spin-axis 
orientation, i.e. PA sp in = PAo. Note that the calculation of 
PA from the Stokes Q and U can only provide a 'headless' 
polarisation vector, P — L exp[tan _1 ([//Q)/(2i)], introduc- 
ing therefore a 180° ambiguity. In addition, as was shown by 
Backer, Rankin & Campbell (1975) and Manchester, Taylor 
& Huguenin (1975), pulsar emission can occur in two orthog- 
onal modes: i.e. the polarisation plane can either coincide 
with the field line-magnetic axis plane or be perpendicu- 
lar to it. Hence, there is an additional 90° ambiguity in the 
projected direction of the spin axis. And so, given the afore- 
mentioned ambiguities, the angle between the spin axis and 
the velocity vector, 9, as is measured on the sky plane (refer 
to Fig.[]J), can only be unambiguously given in the interval 
0°-45°. 

Observationally, there are a number of obstacles that 
prevent us from recovering the intrinsic distribution of spin- 
velocity angles. A major limitation is the unknown veloci- 
ties of the progenitor systems (whether isolated or binary), 
which, if high, can smear even a perfect correlation to an 
apparent misalignment. Very importantly, any study of the 
subject has to rely on the projected velocities and spin-axes 
directions; without knowledge of the pulsar's orientation and 
motion along the LOS, any evidence for alignment may be 
coincidental. There are indeed cases where the inclination 
of pulsar-wind tori can reveal the 3D orientation of the pul- 
sar's spin (e.g. in the case of the Crab and Vela Pulsar Wind 
Nebulae — PWNe; Ng & Romani 2004), but even in those 
cases the unknown progenitor velocity can dominate over 
the measurement errors of the proper motion, depending on 
the distance to the SN site, thus leading to significant uncer- 
tainties on the alignment angles (for a discussion, see §3.2 
of Kaplan et al. 2008). 

2.1.1 The reliability of RVM fits 

There are several examples of pulsars whose profile charac- 
teristics allow one to check the consistency of the derived 
value of PAo from an RVM fit with that based on e.g. pulse 
symmetry, evolution with frequency and features in circular 
polarisation. A few clear examples can be mentioned from 
the work of Johnston et al. (2005, 2007); the reader can 
refer to Fig. [71 for the polarisation profiles of the pulsars 
mentioned here. Success cases of <j>o determination through 
RVM fitting include that of PSR J1740+1311, which has a 
five-component profile with the RVM inflexion point coin- 
ciding with the profile's centre; similarly, PSR J2048— 1616 



has a well-defined RVM swing across its triple-component 
profile, at all examined frequencies, with the inflexion point 
coinciding with mid-point between the outrider components. 

However, there are many cases for which the RVM 
model does not provide a clear determination of <f>o, often 
due to orthogonal-mode emission and other features that de- 
viate from an expected RVM swing; in those cases, e.g. PSR 
J1735— 0724, the profile characteristics may provide alterna- 
tive hints at the location of the minimum approach to the 
magnetic pole. The aforementioned pulsar exhibits a clear 
peak at the location of the central component, where <f>o is 
chosen, but RVM fitting results in large uncertainties due to 
the presence of a hump in the PA profile (most likely caused 
by a transition, a 'jump', between normal- and orthogonal- 
mode emission); in several cases, the inclusion of orthogo- 
nal jumps at preselected locations in the RVM model helps 
towards the determination. Also, it was mentioned earlier 
that the change of handedness of the circular polarisation 
is telling of emission from near the magnetic pole. A vivid 
example is that of PSR J1900-2600, for which the RVM in- 
flexion, the profile's mid-point and the change of handedness 
of the circular polarisation all coincide; thus, for this pulsar, 
the confidence in the determined value of PAo is amongst 
the highest presented here. 

Interestingly, there are cases for which the RVM model, 
even with the inclusion of orthogonal jumps, yields a differ- 
ent solution to that suggested by profile symmetry, etc. A 
good example is that of PSR J0738— 4042, for which a very 
good RVM solution for 0o trails the location of the profile's 
mid-point by fa 4° ; also, a lag of 2° between the profile's mid- 
point and the well-determined RVM inflexion is observed for 
PSR J0452— 1759. Lags of that magnitude may be associ- 
ated with A/R effects, although there are cases where this 
association is tentative (Gupta & Gangadhara 2003). 

Finally, it is also worth noting those pulsars whose PA 
profiles are completely featureless (i.e. flat) across the pulse; 
for such cases, like that of PSR J1915+1009, which has a 
single component, RVM fitting does not help, and profile 
features like the location of the peak intensity are normally 
employed. However, it should be cautioned that such deter- 
mination assumes that the observable emission corresponds 
to a core component, i.e. emanating from near the magnetic 
pole. If this is not the case, the determined cf>o and the cor- 
responding PAo can be completely wrong. 

2.2 Previous Work 

Amongst the recent noteworthy observational investigations 
of pulsar spin-velocity is that of Johnston et al. (2005; here- 
after J05) who used 25 carefully calibrated pulsars, 10 of 
which showed strong evidence for a correlation between spin 
and velocity. The uniform distribution of ^ between 0° and 
45° (i.e. allowing for orthogonal emission) was rejected un- 
der a Kolmogorov-Smirnov (KS) test at the 98% confidence 
level (CL). It should be mentioned that the expected form 
of the *I> distribution for PAo _L PA„ (see Section 14. 1[) was 
also rejected at the 99% CL. In addition, for those pulsars 
in that sample with spin-down ages, r c = P/(2P) < 3 Myr, 
where P is the period time-derivative, the authors found 
that ^ < 10°, although it should be stated that the sample 
was only composed of 4 pulsars. 

The follow up work of Johnston et al. (2007) used a 
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different sample of 22 pulsars, only 7 of which were found to 
be plausibly aligned. The conclusion was that the evidence 
for alignment found in the first study is not supported by the 
separate, follow-up sample. However, it was stressed that the 
pulsars in the second study were both older and more distant 
than those of the first and, therefore, could be more prone 
to the effect of the gravitational potential on the velocity 
direction. 

In an attempt to assess the evidence for spin-velocity 
alignment in view of better statistics, Rankin (2007) re- 
examined the sample of Johnston et al. (2005) and revised 
the PAo values for 23 pulsars of the original sample by 
studying the individual profiles' polarisation and, in particu- 
lar, their frequency evolution. Additionally, Rankin included 
polarisation measurements from Karastergiou & Johnston 
(2006), Xilouris et al. (1991) and Morris et al. (1979, 1981). 
Finally, a small number of additional measurements were 
added to the above sample from new L-band observations 
with Arecibo, amounting to a total of 46 values of PAo and 
the corresponding PA„ values from Hobbs et al. (2005). In- 
terestingly, in this work, the offset angle PAo — PA„ was 
not constrained between 0° and 90° or 0° and 45°, thus 
not folding the polarisation-measurement or orthogonal- 
emission ambiguities into the final result. Nevertheless, the 
distribution of PAo — PA„ (between and 180° ) of the total 
pulsar sample showed an evident clustering of values around 
0° and 90°. Hence, it was argued, based on the majority of 
cases being nearly aligned or nearly orthogonal, that there is 
convincing evidence for a correlation between the spin and 
velocity directions of pulsars. 

Crucially, despite the preferential concentration of val- 
ues around 0° and 90° , in the above studies it was recognised 
that the results are inconclusive as to what SN-kick mech- 
anism is in operation or whether, indeed, the distinction 
between orthogonal or aligned configurations can be trusted 
due to orthogonal- mode emission. 

Finally, even though the data revision by Rankin (2007) 
did not contradict, but rather strengthened, the conclusions 
of Johnston et al. (2005), it highlighted the subjectivity of 
determining PAo through studies of pulsar PA profiles. It 
is therefore fair to say that, ultimately, without the deeper 
understanding of pulsar magnetospheric physics, the con- 
clusions from any such analysis will unavoidably be subject 
to differing interpretations. This is clearly the weak point 
of all polarisation-based analyses to date. Recognising that 
our conclusions may potentially be governed by a bias in the 
selection of PAo, in Section [4.21 we explore the dependence 
of the distribution of "J on various degrees of deviation from 
the published values of PAo . 

2.3 Data Sample 

For this paper, we have tried to compile a database of all 
previously published measurements of PAo that are accom- 
panied by known measurements of PA„ from the literature. 
The aim of this paper is to provide a conclusive statement 
on whether the observed, 2-dimensional spin and velocity 
orientations of pulsars are correlated, for a sample of pul- 
sars with reliable, absolutely-calibrated polarisation mea- 
surements. Although we considered all published measure- 
ments of PAo and PA„, we only included in our analysis 
those measurements for which published polarisation pro- 



files were available and whose quality could be examined. 
The final sample that was used in this work contained 51 
pulsars, whose polarisation profiles have allowed the deter- 
mination of spin-axis directions, and an additional 3 pulsars, 
whose spin-axis directions have been determined through fit- 
ting the the tori of the PWNe seen in Chandra X-ray images 
(see Table [1] for references). 

The pulsar polarisation profiles used in this work came 
from observations that were performed with the Lovell and 
Effelsberg radio telescopes as well as from older observations 
with the Parkes telescope. The details concerning the Parkes 
observations and data selection can be found in the pub- 
lished articles by Johnston et al. (2005, 2007) and Wang et 
al. (2006) . The Lovell observations were carried out in three 
sessions: 2005, 18 September; 2006, 23 March; and 2006, 30 
May. The centre frequency of the observations was 1.4 GHz; 
the 21-cm cryogenic receiver installed at the focus of the 
Lovell telescope is sensitive to both senses of circular po- 
larisation and has a system equivalent flux density (SEFD) 
of 37 Jy. The data were dedispersed coherently and in real 
time using the known pulsar dispersion measure (DM) and 
the COBRA (Coherent On-line Baseband Receiver for As- 
tronomy) coherent-dedispersion back end (Carr 2007). Our 
observations with this backend used different bandwidths 
between 10 and 30 MHz, although the COBRA hardware 
itself is capable of coherently dedispersing up to 100 MHz. 
Depending on the pulsar luminosity, the integration time of 
those observations varied between 3 and 15 minutes. A total 
of 12 pulsars were observed with the Lovell system. 

The Effelsberg observations were performed at two fre- 
quencies: 2.7 GHz and 4.85 GHz, utilising 80 and 500 MHz of 
bandwidth, respectively. The 2.7-GHz receiver is equipped 
with circular-polarisation feeds and provides an SEFD of 
11 Jy; the 4.85-GHz receiver also has circular-polarisation 
feeds and an SEFD of 18 Jy. The observations took place 
during 2006, 12-14 March and 2006, 27 May, during which 
a total of 32 pulsars were observed, with 21 of these be- 
ing different from those observed with the Lovell. The data 
were acquired and folded with the Effelsberg Pulsar Ob- 
servation System (EPOS) and the PUB86 pulsar back end 
(Kramer 1995). Dedispersion was not necessary, as the dis- 
persive delay for the observed pulsars, given the available 
bandwidth at the relevant observing frequencies, was neg- 
ligible compared the pulse period: i.e. the dispersive delay 
between two frequencies, /i and /2 (in MHz), for a pul- 
sar with a dispersion measure DM (in cm~ 3 pc), is given by 
At = (4.15 x 10 6 ms) x DM x (/f 2 - f^ 2 ); for a bandwidth 
of 30 MHz at 2.7 GHz and at 4.85 GHz, this delay is « 12 
/xs and « 2 fis per DM unit, respectively. The pulsars in 
our sample have DM < 100 cm _3 pc and P > 0.2 s, which 
means that dispersion smearing is confined to, at most, only 
~ 0.5% of the pulse period. 

In both the Lovell and Effelsberg polarisation obser- 
vations, the Stokes data of each observation were abso- 
lutely calibrated by correcting, amongst other effects, for 
the parallactic-angle rotation during the observation and 
the Faraday rotation, using the published rotation measure 
(RM) of the pulsar (Manchester et al. 2005). In general, the 
error on RM, ctrm, could significantly impact on the PAo 
error, <tpa; the latter is given by ctpa = c 2 ctrm// 2 , where c 
is the speed of light and /, the observing frequency. In John- 
ston et al. (2005, 2007), the RM errors were folded into the 
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Table 1. Table of velocity and spin orientation for the 54 pulsars used in our study: Column 3 shows the position angle of the projected 
velocity vector, PA„, derived from pulsar timing and VLBI (Very Long Baseline Intcrferometry) observations; Column 5 lists the position 
angle of the spin-axis projection onto the sky plane, PAo, as was derived from radio polarisation measurements and fits to the PWN 
X-ray images; finally, Column 7 shows the minimum offset angle between projected velocity and spin, \E" = (PAo — PA„) mod 45°, which 
accounts for polarisation and orthogonal- mode ambiguities. The statistical errors on the PA V , PAo and are shown in parentheses 
and refer to the last significant digit. More information on the observations and the respective techniques used can be found in the 
publications referenced in Columns 4 and 6. 



PSR 


logTc 


PA„ 


Reference 


PAo 


Reference 






[yr] 


[°] 




[°] 




[°] 


J0139+5814 


5.6 


30(13) 


[1] 


-71(5) 


[18] 


11(14) 


J0152-1637 


7.0 


173(3) 


[5] 


92(9) 


[19] 


8(9) 


J0304+1932 


7.2 


-9(11) 


[2] 


12(10) 


[18] 


21(15) 


J0332+5434 


6.7 


-61(1) 


[3] 


-4(5) 


[18] 


33(5) 


J0358+5413 


5.8 


69(16) 




-33(5) 


[18] 


12(17) 


J0452-1759 


6.2 


72(23) 


[5] 


47(3) 


[20] 


25(23) 


J0454+5543 


6.4 


-72(3) 


[1] 


-66(11) 


[18] 


6(11) 


J0525+1115 


7.9 


132(16) 


[5] 


-65(4) 


[21] 


17(16) 


J0534+2200 


3.1 


292(10) 


[6] 


124.0(l) b 


[22] 


12(10) 


J0538+2817 


4.6 


-24.1(2) 


[7] 


-49(10) 


[18] 


25(10) 


J0543+2329 


5.4 


57(20) 


[1] 


-100(5) 


[18] 


23(20) 


J0630-2834 


6.4 


294(3) 


[8] 


26(2) 


[21] 


2(4) 


J0659+1414 


5.0 


93.1(4) 


[9] 


-86(2) 


[20] 


1(2) 


J0738-4042 


6.6 


313(5) 


[8] 


-21(2) 


[20] 


26(5) 


J0742-2822 


5.2 


278(5) 




-81.7(1) 


[21] 


0(5) 


J0820-1350 


7.0 


159(6) 


[5] 


65(2) 


[21] 


4(6) 


J0826+2637 


6.7 


-34(1) 


[2] 


36(10) 


[18] 


20(10) 


J0835-4510 


4.1 


301.0(1) 


[10] 


36.8(1) 


[21] 


5.8(1) 


J0837+0610 


6.5 


2(5) 


[5] 


18(5) 


[20] 


16(7) 


J0837-4135 


6.5 


187(6) 


[8] 


-84(5) 


[20] 


1(8) 


J0922+0638 


5.7 


12.0(1) 


[11] 


-56(5) 


[18] 


22(5) 


J0953+0755 


7.2 


355.9(2) 


[3 

L -'J 


14.9(1) 


[21] 


19.0(2) 


J1136+1551 


6.7 


348.6(1) 


[3] 


-78(2) 


[21] 


23(2) 


J1239+2453 


7.4 


295.0(1) 


[3] 


-66(1) 


[21] 


1(1) 


J1430-6623 


6.7 


236(9) 


[12] 


-28.5(7) 


[21] 


6(9) 


J1453-6413 


6.0 


217(3) 


[12] 


-56.9(4) 


[21] 


4(3) 


J1456-6843 


7.6 


252.7(6) 


[13] 


-31.6(6) 


[21] 


14.3(8) 


J1509+5531 


6.4 


47(2) 


[2] 


-3(7) 


[18] 


40(7) 


J1604-4909 


6.7 


268(6) 


[13] 


-17(3) 


[20] 


15(7) 


J1645-0317 


6.5 


353(3) 


[8] 


56(4) 


[21] 


27(5) 


J1709-1640 


6.2 


192(16) 


[5] 


15(2) 


[21] 


3(16) 


J1709-4429 


4.2 


160(10) a 


[14] 


163.6(7) b 


[23] 


4(10) 


J1735-0724 


6.7 


355(3) 


[8] 


55(5) 


[20] 


30(6) 


J1740+1311 


6.9 


228(6) 


[8] 


-46(4) 


[21] 


4(7) 


J1801-2451 


4.2 


270 c 


[15] 


-55(5) 


[20] 


35(5) 


J1820-0427 


6.2 


338(17) 


[5] 


42(3) 


[20] 


26(17) 


J1844+1454 


6.5 


36(15) 


[5] 


-52(2) 


[21] 


2(15) 


J1850+1335 


6.6 


237(16) 


[5] 


-45(3) 


[20] 


12(16) 



final, published uncertainty on PAo; this was not done for the 
pulsars of observed with the Lovell and Effeslberg telescopes. 
However, for the vast majority of those pulsars <trm < 1 rad 
m , which at the observing frequencies in question meant 
that o pa < 1°. More importantly, for all those pulsars the 
measurement errors on PA„ and PAo dominate over the er- 
rors due to RM uncertainties: the dominant measurement 
errors are due to (a) the limited signal-to-noise ratio (S/N) 
in linearly polarised intensity, L, i.e. <tpa ~ 07 /L, where 07 
is the rms of the total intensity; and (b) due to uncertainties 



in proper motion determination. The worst-case scenario in 
our sample is PSR J0922+0638, observed with Effelsberg at 
2.7 GHz: this pulsar has RM = 32 ± 2 rad m~ 2 , which cor- 
responds to <tpa ~ 1?5 and which is still small compared to 
the combined measurement error on ^, i.e. <r* =5°. For a 
further discussion on the systematic effect of erroneous de- 
terminations of PAo in statistical studies of pulsar spinave- 
locity alignment, the reader can refer to Section f4. 21 

Finally, the frequency-independent, systematic PA ro- 
tations, caused by phase offsets during the processing chain, 
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Table 1. Continued. 



PSR 


logr c 


PA„ 


Reference 


PA 


Reference 






[yr] 


[°] 




[°] 




n 


J1900-2600 


7.7 


202.8(7) 


[16] 


-43(2) 


[21] 


24(2) 


J1913-0440 


6.5 


166(11) 


[5] 


-68(2) 


[21] 


36(11) 


J1915+1009 


5.6 


174(15) 


[5] 


85(3) 


[20] 


1(15) 


J1921+2153 


7.2 


34(12) 


[5] 


-35(2) 


[21] 


21(12) 


J1932+1059 


6.5 


65.2(2) 


[11] 


-11.3(1) 


[21] 


13.5(2) 




u.u 


i-IO(l) 


[51 


1 1 (7) 


[911 


i4n ~i 


J1937+2544 


6.7 


220(9) 


[5] 


-9(5) 
85(5) b 


[20] 


41(10) 


J1952+3252 


5.0 


252(7) 


[17] 


[22] 


13(9) 


J1955+5059 


6.8 


-27(3) 


[5] 


19(10) 


[18] 


44(10) 


J2013+3845 


5.6 


47(11) 


[5] 


83(5) 


[18] 


36(12) 


J2018+2839 


7.8 


23(2) 


[3] 


89(20) 


[18] 


24(20) 


J2022+2854 


6.5 


11(1) 


[3] 


12(10) 


[18] 


1(10) 


J2022+5154 


6.4 


-24.5(9) 


[3] 


13(5) 


[18] 


38(5) 


J2048-1616 


6.5 


92(2) 


[4] 


-13(5) 


[20] 


15(5) 


J2157+4017 


6.8 


-81(3) 


[1] 


79(10) 


[18] 


21(10) 


J2330-2005 


6.7 


86(2) 


[8] 


21(10) 


[20] 


25(10) 



References — Proper motions: [1] Harrison et al. (1993); [2] Lyne et al. (1982); [3] Brisken et al. (2002); [4] Fomalont et al. (1997); 
[5] Hobbs et al. (2005); [6] Caraveo & Mignani (1999); [7] Ng et al. (2007); [8] Brisken et al. (2003a); [9] Brisken et al. (2003b); 
[10] Dodson et al. (2003); [11] Chatterjee et al. (2001); [12] Bailes et al. (1990b); [13] Bailes et al. (1990a); [14] Wang et al. (2006); 
[15] Blazek et al. (2006); [16] Formalont et al. (1999); [17] Migliazzo et al. (2002). Polarisation: [18] Carr (2007); [19] Wang et 
al. (2006); [20] Johnston et al. (2007); [21] Johnston et al. (2005); [22] Ng & Romani (2004); [23] Romani et al. (2005). 

a. Assumes an association between the pulsar and the SNR G343. 1—2.3 (see e.g. Dodson & Golap 2002). 

b. Spin axis direction derived from fitting the tori of the PWN seen in Chandra images. 

c. Position angle corresponding to the upper limit on the westerly motion of this pulsar, derived from interferometric observations. 



were removed by comparing the PAs of absolutely calibrated 
Parkes observations of 5 bright, polarised pulsars (covering 
a wide range of parallactic angles) to those obtained with 
the Lovell, after the aforementioned calibration steps. Sim- 
ilarly, for the Effelsberg observations, PSR B 1929+ 10 was 
used as a polarimetric calibrator. Further details concerning 
the new observations can be found in Carr (2007). 

The statistical investigation of the relative phases of 
the pulsar spin and velocity vectors is strongly affected by 
the systematic uncertainties (see Section f 4 . 2 [I . Therefore, a 
careful selection of the best constrained PAs is necessary 
to minimise the ambiguity of our conclusions. One major 
systematic uncertainty comes from the determination of the 
spin axis direction, which is most commonly derived from 
pulsar polarisation profiles. Fig. [7] shows all the integrated 
polarisation profiles, from the aforementioned observations, 
that were used in this work. The lower panel of each plot 
shows the total intensity profile (black) and the linear (red) 
and circular (blue) polarisation profiles. The top panels show 
the profile of the polarisation position angle (PA) across the 
pulse. In each of the plots a vertical, dashed black line is 
drawn at PAo, i.e. the position angle of polarisation at the 
point of closest approach of the observer's LOS to the mag- 
netic pole. 

As well as the polarisation measurements, Fig. [7] shows 
in green horizontal dashed lines, the position angle of the 
proper motion, PA„, in order to provide a direct compari- 
son between the orientations of the spin and velocity vectors. 
The components of pulsar proper motions used to calculate 



PA„ are the most accurate measurements available, taken 
from Hobbs et al (2005), Johnston et al. (2005) and other 
sources referenced in Table Q] Given the measured values 
of PA„ and PAo, we can calculate the offset angle, "J, be- 
tween the sky-projected spin-axis orientation and velocity 
direction. Taking into account all the ambiguities mentioned 
above, ^ is given by $ = (PAq — PA„) mod 45°, with astan- 
dard deviation calculated from cr* = \J <Tq + a% , where oo 
and <r„ are the standard deviations of PAo and PA„, respec- 
tively. Here and later in our analysis, we assume Gaussian 
statistics for the PA distributions. All PA„, PAo values and 
the corresponding $ angles for the pulsars in our sample are 
shown in Table [T] 

In summary, we have tried to address the open ques- 
tion of whether pulsar spins and their respective velocities 
are preferentially aligned, using 54, good-quality (PAo, PA„) 
value pairs, of which 16 came from new measurements by 
Carr (2007). To date, this is the largest sample of pulsar 
spin-velocity offset angles considered in such investigation. 



3 DATA ANALYSIS 

If there is a SN-kick mechanism that promotes the align- 
ment between the pulsar spin and the kick velocity, then for 
a sample of pulsars that possess this property one should ex- 
pect to see the majority of ^ values distributed towards 0°. 
On the other hand, a uniform distribution of offset angles 
would indicate no correlation. A complication comes from 
the fact that the gravitational potential of the Galaxy alters 
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Figure 2. (a) Distribution of vp, for the pulsars of Tablc[T] from 10 e Monte Carlo simulations of PAo and PA„, based on their published 
error values. The vertical, dotted grey line marks the mid-range of the parameter space = 22?. 5). (b) Distribution of the KS statistic, 
D, from the above simulations, from the comparison of the distribution shown in (a) with a uniform distribution. The probability of 
rejecting the uniform distribution, p(D), and the lcr confidence interval is also shown, (c) Two-dimensional PDF of PA„ and PAo, 
generated from random Gaussian iterations of the values in our data within their respective errors. To produce this map, we re-defined 
PA V as PA„ — PAq mod 90° , considering only the polarisation ambiguity, and forced PAo to 0. The final distribution shown was generated 
by randomising the the re-defined PA V and PAo within their errors and wrapping their difference modulo 90° , when necessary. The solid 
grey and dashed colour lines overlaid on the map correspond to the 0°,±30°,±60° and ±90° spin-velocity offset contours. The colour 
scale indicates the probability, p, increasing from blue towards red with graduated colours, corresponding to the different points on the 
map. 



the pulsar velocity vectors during the lifetime of a pulsar. 
Hence, such correlation would be stronger for young pulsars 
(e.g. < 1, 000 kyr) and vanish for older pulsars (e.g. > 10, 000 
kyr) . The investigation of spin-velocity alignment as a func- 
tion of age is an interesting problem in its own right, but it 
critically depends on the reliability of pulsar ages and dis- 
tances. In the present paper, we do not make a distinction 
between different pulsar age groups or consider the effect of 
the Galactic gravitational potential on pulsar velocities; a 



detailed study of the above problem will be presented in a 
separate, follow-up publication. 

A simple way of checking whether the conclusions de- 
rived from the above studies hold for our selected sample 
of pulsars is by quantifying how different the distribution of 
$ is from a uniform distribution, using a statistical test for 
similarity, like the KS test. Our limited sample of N = 54 
values of *P can be simply binned in steps of 5°, between 
and 45° , to produce a discrete distribution of the offset 
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angles. However, as is clear from Table [TJ the errors on PAo 
and PA„ are significant and have to be taken into account in 
quantifying the uniformity of the distribution. Statistically 
speaking, a parent population of pulsars, part of which is 
our sample, would produce a distribution of ^ that is the 
convolution of the individual Gaussian distributions of each 
pulsar in our sample, each having a standard deviation equal 
to the measured one. A numerical way of approximating the 
parent distribution using the available data is by perform- 
ing a large number of Monte Carlo (MC) iterations, allowing 
in each random instance the values of PAo and PA„ to be 
distributed normally in A r (PAo, erg) and N(PA V , a v ), respec- 
tively. The resulting values of ^ from each iteration are then 
binned in \I/ £ [0°,45°]. Following this recipe, we generated 
10 6 histograms of ty, which we then averaged to produce the 
final, aggregate distribution of the distribution is shown 
in Fig. [2^,. This MC procedure shuffles the data positions be- 
tween bins of ^, according to the er* of each pulsar, which in 
turn causes the fluctuation of the bin heights, rik, between 
iterations. Therefore, a good measure of the statistical error 
on the feth bin, <7fc, in the aggregate distribution is the stan- 
dard deviation of the bin's height across all MC iterations: 

i-e. (Jfc = V( n fe)MC - (rifc)MC- 

We examined the uniformity of the aggregate distribu- 
tion of ^ that was generated from the sample of Table Q] us- 
ing the described procedure. At first glance, the distribution 
appears non-uniform with ~ 70% of the ^ values distributed 
below the median (= 22? 5). As mentioned earlier, a suitable 
test for quantifying the similarity of the aggregate distribu- 
tion to a uniform one is the KS test. In each of the iterations 
of the MC procedure, we performed a KS test between the 
randomly generated distribution of *]/ and a uniform prob- 
ability function, p(^) = 1/45 deg -1 . Thus, we generated 
10 6 values of the KS statistic, D. As is shown in Fig. 03, 
the values of the statistic D are distributed in a Gaussian 
fashion, which justifies choosing the mean, D, and the stan- 
dard deviation, an, as the most probable value and its 68% 
CLs, respectively. Based on these values, the probability of 
rejecting uniformity under the KS test was calculated from 
D, i.e. p(D), and the probabilities corresponding to the 68% 
CLs, from D±o~d, i.e. p(D^fao jE These values are shown in 
the legend of Fig.03. Hence, the apparent non-uniformity of 
the $ distribution is confirmed by the KS test, which rejects 
uniformity at the 99.87_J'2% level, i.e. at roughly 3<r. 

In addition to plotting the histogram of the wrapped 
angle ^, we wished to examine the raw, 2-dimensional dis- 
tribution of PAs, i.e. PA„ and PAo, without considering the 
possibility of orthogonal emission but still allowing only the 
polarisation-measurement ambiguity. Since we are only in- 
terested in the relative direction of the spin and velocity 
axes, we set all measured PAo values to 0° and re-calculated 
PA„ as PA„ — PAo, which was forced to be between —90° 
and 90°. By randomising, as before, both angles within their 
respective statistical errors (and wrapping their difference 
modulo 90°, as necessary), we generated a large number 
of (PAo,PA„) pairs. The 2-dimensional probability density 
function (PDF) of the generated distribution is shown in 
Fig- [It- The PDF clearly shows an appreciable concentration 



Note the inverted signs, as we are interested in the probability 
of rejecting uniformity. 



of values around (0°, 0°) and near (0°, ±90°): this picture is 
consistent with the 1-dimensional distribution of ^, as all 
the values inside the most significant contours (p > 0.0015) 
would fall into the bins below the mid-range of >3/, if folded 
back into the 0° — 45° interval. Also, at first glance, the cen- 
tral island near 0° appears the most significant of the three, 
having a central probability contour (p > 0.002) that is more 
extended than both the one near —90° and that near 90°. 
However, it is difficult to visually compare the likelihood of 
the aligned (~ 0°) with that of the orthogonal (~ ±90°) con- 
figurations, only based on the data, without modelling their 
expected distributions. A quantitative comparison between 
these two configurations, and for various PA^ distributions, 
is presented in the following Section. 

In conclusion, if we take the PDF of Fig.[5Ji at face value, 
it appears that the concentration of ^ values below 22? 5 in 
the histogram of Fig. [2^, is the result of not only nearly par- 
allel but also nearly orthogonal spin-velocity configurations. 
In particular, if we consider that having PA„ — PAo = 90° 
and —90° is essentially the same configuration, then both 
aligned and orthogonal configurations appear to have com- 
parable probability of occurrence, based on the PDF shown. 
Equally interesting is the fact that PA offsets in the ranges 
[—60°,— 30°] and [30°, 60°] appear very unlikely according 
to the PDF, which confirms that the observed alignment is 
robust. However, as was conclusively shown with radio and 
X-ray observations of the Vela pulsar (Ng & Romani 2004; 
Ng & Romani 2007) , configurations that appear orthogonal, 
due to orthogonal-mode emission, may actually be aligned. 
In such case, it is possible that Fig. [2t shows a distribu- 
tion of only aligned configurations, some of which appear 
orthogonal; while of course the opposite is also possible. 



4 DISCUSSION 

4.1 Toy- Model Simulations 

4-1.1 Model description & geometry 

In order to further quantify the margins of the process that 
could potentially lead to the observed non-uniformity in the 
distribution of we performed numerical, toy-model simu- 
lations. We simulated a high number of pulsar-velocity direc- 
tions over 47t rad 2 , assuming a parent process that results in 
(a) a peri-axial distribution of velocities (i.e. case of aligned 
kicks), (b) a distribution of velocities nearly aligned with 
the equatorial plane (i.e. orthogonal kicks) and (c) a uni- 
form distribution of velocities with no preferred direction. 
Furthermore, in order to simulate the random orientations 
of pulsar spins with respect to the observer, we allowed the 
observer's LOS to be distributed uniformly over the entire 
sky but kept the pulsar spin-axis orientation fixed. After pro- 
jecting the spin and velocity directions onto the observer's 
plane (perpendicular to the LOS), we folded the resulting 
values of ^ in the 0° — 45° interval (as was done with the 
real data) . The character of the resulting distributions of \f r 
for those three cases can be seen in Fig. [3] where we arbi- 
trarily used a Gaussian standard deviation of 10° for the 
PA„ distributions (hereafter simply referred to as a v ). 

The observed distribution of *t can be more clearly un- 
derstood with the help of the schematics of Fig. [4] which 
show the three-dimensional geometry invoked in our simu- 
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Figure 3. Simulated distributions of \1/ = (PAo — PAi,) mod 45°, for three cases of alignment between the pulsar velocity, v, and its 
spin axis, S: (a) the distribution of v is aligned with S to within a v \ (b) the distribution of v is orthogonal to S to within tr„; and 
(c) the distribution of v is uniform on a sphere. The geometric details of the simulation are explained under Fig. [4] For case b, the ^ 
distribution for cr v = 0° is shown with a dashed white line. 



lation. It is clear that a preferred alignment would lead to 
the ^ distribution of Fig. [3^, where > 90% of the values 
are below 20° , and with the distribution practically vanish- 
ing towards = 45° . The underlying geometry of this case 
is shown in Fig. [4^,, configuration (i). It is also worth not- 
ing that uniformly distributing the observer's LOS over 47t 
does not significantly affect the shape of this distribution, 
because the only possibility of having ^> 0°, in this case, 
would be for LOS that are nearly parallel to S, which is a 
configuration with low probability. Such a low-probability 
case is shown in Fig. |4p. 



On the other hand, the orthogonal case results again in 
an uneven distribution with an excess of values below 20° . As 
can be seen in Fig.|4^, configuration (ii), the projection of the 
three-dimensional velocity distribution onto the observer's 
plane, for an arbitrary angle between K and S, is an ellipse 
with its minor axis aligned with S± (i.e. the projection of S 
onto the observer's plane). If we allow the distribution of v to 
have a width around the equator, then the projection of this 
equatorial band onto the observer's plane is an elliptical ring- 
shaped area — formed by the superposition of ellipses within 
the equatorial band — that is broadest along the direction 
of S± and whose width diminishes perpendicular to it. At 
the extreme ends of the possible — KZS configurations, 
that projection degenerates into either a diametral band (for 
6 — 90°; see Fig. 2t) or a circular band (for 8 = 0° or 
180°; not shown). In the former case, the most probable 
values of $ are ~ 90° , whereas in the latter case \t can take 
all values in [0°,180°] with equal probability. Crucially, as 
was mentioned in the previous paragraph, the probability of 
having 9 = 0° (or 180°) is significantly lower than that for 
6 = 90° , in a uniform distribution of LOS and S directions. 
Therefore, after folding all angles modulo 45°, there will 
be an excess of angles near *t ~ 0° compared to those near 
45° . And this is true independently of any dispersion of the 
velocity directions around the equator, as is clearly shown 
with the dashed white line in Fig. [3]d, where no deviation of 
the velocity direction from orthogonality was assumed. 



Finally, as expected, the case of uniformly distributed 
velocities results in a uniform distribution of \&. 



4-1.2 Fits to the data 

The significant departure from a flat distribution of the ob- 
served distribution of \& gives us the opportunity to investi- 
gate the properties of the parent distribution of 3D velocities 
by comparing it with our toy-model simulation. In order to 
produce a quantifiable measure of similarity, we performed 
a two-sample KS test between the observed and simulated 
distributions of ty, for a range of a v values. In addition to the 
free parameter a v , we added a fixed amount of observational 
error on ^, roughly equal to the average value measured in 
our data set: i.e. <r* = 10°. Thereafter, we generated 100 
simulated distributions of ^ for a v € [0°, 45°], each contain- 
ing 100,000 values. For each of the simulated distributions, 
we generated 10,000 MC data sets by normally distributing 
the measured $ values according to their standard devia- 
tions, similarly to the procedure done for Fig. [2] 

For each a v , we performed a two-sample KS test be- 
tween the simulated distribution ^(a v ) and each of the 
10,000 data sets. Each of these tests produces a KS statis- 
tic, D\ the whole procedure produces a distribution of 10,000 
values of D, similar to that of Fig. [3> — which was derived 
from the one-sample KS test between p(^) = 1/45 deg -1 
and a large number of data sets. From that distribution, 
we calculated the mean KS statistic, D, and its 68% CLs, 
D ± od- Finally, the corresponding probability of rejecting 
the null hypothesis under the KS test was calculated from 
D and D ± an- i.e. p(D; a v \Ho) and p(D ± o\d; a v \Ho), re- 
spectively. The resulting PDFs for the cases of aligned and 
orthogonal kicks are shown in Fig. As well as the central 
PDFs, which correspond to the p(D; <r v \Ho) values (black 
lines), we have also plotted those corresponding to the 68% 
CLs on D (dashed-dotted lines), as was explained above. For 
clarity, the areas between p = and each of the PDF func- 
tions shown (with solid and dashed-dotted lines) have been 
shaded with dark-red and dark-blue colours, for the central 
PDFs of the aligned and orthogonal case, respectively, and 
with light-red and light-blue colours, for the PDFs of the 
aligned and orthogonal case, respectively, corresponding to 
the 68% limits on D. An immediately evident feature of the 
generated PDFs is that they show clear global maxima for 
certain values of cr„. For the case of alignment, the central 
PDF gives a 68% confidence interval around the maximum 
of a v — 26lg deg. Interestingly, the PDF in the orthogonal 
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Figure 4. Schematics of the geometry used in our spin-velocity simulation code: (a) the top schematic shows the general case, where 
the observer's LOS, K , and the spin-axis vector, S, form an arbitrary angle, 9 = KZS; (b) the middle schematic shows the special case, 
where K/LS is small, i.e. 9 <C 45°; (c) finally, the bottom schematic shows another special case, where 9 = 90°. Assuming a uniform 
distribution of the observer's LOS and the spin-axis directions over 47t, case (b) is less probable than case (c). In all schematics, the real, 
three-dimensional vectors for the pulsar spin and a few exemplifying velocity vectors, v, are represented with solid, black vectors; the 
solid, grey vectors correspond to the projections of the three-dimensional vectors onto the plane perpendicular to K (observer's plane). 
Plot (a) showcases two configurations of the pulsar's velocity vectors with respect to K and S, designated with the superscripts (i) 
and (ii): those of (i) are part of a distribution of velocities about the spin axis, within a solid angle surrounding S (probability density 
contours of increasing angular separation from S are also shown); those of (ii) are part of a distribution of velocities along the equator 
(dashed black line), within an equatorial band (delimited by solid black lines). The projection of the three-dimensional distributions of 
velocities onto the observer's plane, for both configurations, is shown as two-dimensional, grey-shaded probability density contours (for 
configuration (ii), the projection of the distribution that is below the observer's plane is also drawn). On the observer's plane, the angles 
'J between the projected spin and velocity vectors, S±_ and v±, respectively, are also shown. Using the same conventions as in (a), plot 

(b) shows the case where the velocities are distributed about the pulsar's spin axis and are also roughly aligned with K. Finally, plot 

(c) shows the special case where the velocities are distributed along the equatorial plane of the pulsar and the pulsar's spin axis lies in 
the observer's plane; in this case, the projection of the velocity distribution degenerates into a diametral band, yielding high probability 
for * ~ 90°. 
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Figure 5. The PDFs of <r„, in the case of aligned (red-shaded areas) and orthogonal (blue-shaded areas) kicks, from the comparison 
(under a two-sample KS test) of the data distribution of for all pulsars in our sample with the simulated histograms of Fig. [3] for a 
range of trial <r v values. The central PDFs, for each case, are delimited with solid black lines and correspond to the mean value of the KS 
statistic, D, from our MC simulation, for each <j v : i.e. p(D). The PDFs delineated with dash-dotted black lines correspond to the 68% 
confidence limits of the KS statistic, D ± 073: i.e. p(D ± crp). The correspondingly coloured error bars along the horizontal axis denote 
the median and lcr asymmetric errors on a v , based on the central PDFs of the alignment and orthogonal case. Finally, the black error 
bar corresponds to the step in cr v used in the simulated \l> distributions. 



case also exhibits a clear peak, but in that case the peak 
probability is less than half that of the case of alignment 
(p± « 0.25lo.?s); the latter case results in a probability of 
60 ± 30% that the observed and simulated distributions are 
the same. For the central PDF, the case of orthogonal kicks 
gives a peak at a v — 3^2 deg. 

4-1-3 Conclusions 

In conclusion, the comparison between our toy-model sim- 
ulations and the data have revealed a preference towards 
alignment between velocity and spin, if a v ~ 30° — al- 
though it should be stressed that in both cases the PDFs 
are very broad and irregular, and therefore a v is not very 
well constrained. This result also reveals that, in fact, the 
central island of Fig. [2}: is more significant than those close 
to ±90° . The most probable ranges of o~ v obtained from our 
analysis can be compared with those derived from the sim- 
ulations of KPP using data from Ng & Romani (2007) and 
Rankin (2007; see Introduction). For the case of alignment, 
our range of a v values implies a much larger spread of kick 
angles than the narrow distribution of kick angles using the 
data by Rankin (< 5° — 10°). On the other hand, the case 
of orthogonal kicks, which is favoured by scenarios of binary 
disruption, is much more consistent with the above, narrow 
range. Interestingly, according to KPP, the data by Rankin 
can hardly be fitted without including binary progenitors in 
the kick model. Furthermore, as was mentioned in the in- 
troduction, the data of Ng & Romani can be fitted by a 



much broader distribution of kick angles (< 20°). This re- 
sult is more consistent with our case of alignment, although 
we find a somewhat broader spread of angles (< 20° — 35°). 

Ultimately, our analysis shows that both aligned and 
orthogonal kicks appear probable and that it is difficult to 
confidently rule out one of the modes in favour of the other. 
Hence, we are forced to consider the possibility that an en- 
semble of kick mechanisms may be in operation in the SN 
population, resulting in a mixture of aligned and orthogonal 
kicks. Also, the significant fraction of pulsars that are born 
in binary systems complicates things further: our analysis 
cannot discriminate between spin-velocity angles that are 
the result of a binary motion superimposed on the SN birth 
kick. Furthermore, as will be discussed in a follow up paper, 
the role of the Galactic gravitational potential is also an 
important one and can drastically distort any spin-velocity 
correlation that may have been originally in place. 

As a final remark, it should be stated that our simu- 
lation has not considered the pulsar-beam orientation with 
respect to the observer, i.e. the detectability of the simu- 
lated pulsars. In order to investigate whether beam geome- 
try has an effect on the simulated distributions, we reran our 
simulation, this time excluding pulsars whose beam did not 
intersect with the observer's LOS. This condition demands 
certain assumptions about the angular radius of the pulsar 
emission cone, p, and the distribution of the magnetic incli- 
nation, a. Based on the work of Lyne & Manchester (1988), 
we used p = 6?5P~ 1//3 , which corresponds to a median of 
p ~ 5°, for our pulsar sample; hence, we chose p = 5° as 
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the beam width. Furthermore, there are a number of PDFs 
for a, p(a), in the literature (see e.g. Gil & Han 1996): the 
most prominent ones are a uniform distribution in a and 
in cos a, i.e. p(a) — 2/tt and p(ot) = sin a, respectively — 
where a £ [0, 7t/2]; the latter of the two distributions is 
just a uniform distribution of the beam directions over 4n 
(cf. the same distribution was assumed for K). Other, more 
complicated distributions exist, but, in general, the limited 
sample of available measurements combined with the afore- 
mentioned issues with RVM fitting (and other methods of 
obtaining a; e.g. Rankin 1990), as well as our lack of un- 
derstanding of the evolution of a with pulsar age (e.g. Tau- 
ris & Manchester 1998; Weltevrede & Johnston 2008), does 
not warrant in our opinion such complicated descriptions. 
As with the observer's LOS, we randomly chose a from the 
uniform distribution, accepting only configurations for which 
a — p ^ 8 ^ a + p or (7r — a) — p ^ 9 ^ (tt — a)+p. As before, 
a v was kept fixed to 10° . As is expected for a uniform distri- 
bution in a, for each randomly generated configuration of K 
and v, only a certain fraction of beam orientations are ac- 
cepted; but that fraction does not depend on a and, hence, 
the final $ distributions remain unchanged, in this case. On 
the other hand, choosing a uniform distribution in cos a re- 
sulted in practically the same \l/ distribution for the aligned 
case. For the orthogonal case, a uniform cos a distribution 
resulted in a slightly steeper distribution for compared 
to the one shown in Fig. \5jp. However, by increasing a v by 
20%, we were able to match the two distributions, meaning 
that within 20% adding pulsar detectability does not affect 
the simulated distributions. In summary, the magnitude of 
the errors on the determined values of a v , for the orthogonal 
case, ranges from 60% to 260%; this dominates over the bias 
from excluding detectability. And since there is currently 
no consensus in the literature as to which distribution of a 
gives the best description, we decided not to take this into 
account in our simulations. 



4.2 Systematics Inspection 

The above treatment neglects possible systematic errors, 
which can have a significant effect on the distribution of 
Fig. [2] It is therefore important to test the robustness of 
the observed correlation between the spin and velocity posi- 
tion angles against systematic departures from the data set 
of Table [T] The following paragraphs consider two, poten- 
tially significant sources of bias: the choice of PAo and the 
influence of the J05 subset on the degree of correlation. 

4.2.1 The choice of PA 

As was stated in Section ^. 21 a potential bias in our study of 
pulsar-spin alignment using pulsar polarisation is the choice 
of PAo. The determination of the reference PAs that were 
used in this paper was based on a number of criteria em- 
ployed by the respective investigators (see Section [2]). How- 
ever, due to the limited statistics, the results of such anal- 
ysis are potentially sensitive to erroneous determinations of 
PA . For example, the PA profile of PSR J0630-2834 sweeps 
across the entire range of possible angles, i.e. APA PP ~ 180°; 
here, APA PP is defined as the peak-to-peak difference of the 
PA profile for across all phases with L/ai > 3, i.e. S/N in L 



greater that 3 (as shown in Fig. [7|). Clearly, the decision of 
PAo for this pulsar will have a significant impact on whether 
it shows spin-velocity alignment or not; on the other hand, 
the characteristic RVM shape of the PA profile assists an 
easier determination of the inflexion point. In contrast, PSR 
J2013+3845 possesses a flat PA profile (APA PP ^30°), so 
an error on the choice of PAo would not significantly alter 
the result — but whether the relative constant PA across 
the profile corresponds to the PA at the minimum approach 
to the magnetic pole or not is unknown. In order to test how 
sensitive our data sample is to the choice of PAo, we gener- 
ated a high number of data sets, choosing each time a dif- 
ferent PAo for each pulsar, from N(PAo, qAPA pp ), where q 
ranged from 1 (corresponding to a standard deviation equal 
to the entire range of possible PAs) to 1 /w (corresponding 
to a standard deviation equal to 10% of APA PP ). The value 
of APA PP varied from pulsar to pulsar according to the PA 
profiles of Fig. [7] Again, for each of the above data sets we 
performed the usual MC procedure, in order to estimate the 
mean values of the KS statistic. The standard deviations 
of the randomly generated values of PAo, from the previ- 
ous step, were kept equal to those of the published PAs, 
i.e. a = (To, independently of the actual value that a PA 
may have had at those same phases (PAs near the edges of 
the profiles have in general higher a due to the lower S/N). 
The probabilities corresponding to the median value of the 
KS distribution of means and its 68% CLs are shown as a 
function of q in Fig. [6^. 

The most apparent feature of Fig. [EJi is the evident 
trend towards higher probabilities for decreasing errors on 
PAo. It is true, however, that the degree of spin- velocity 
alignment appears very sensitive to the magnitude of the 
departure from the published values of PAo: standard de- 
viations of > 30% of the APA PP result in practically no 
alignment (p < 0.8), whereas rejection of Ho at the 2a level 
is only achieved for q < 0.1. The immediate conclusion from 
this test is that the observed alignment remains robust to 
within a 10% systematic error on the choice of PAo (for a 2a 
rejection of _ffo)- As was suggested above, beyond the impact 
of systematic errors induced by human choice, such tests are 
unable to account for the intrinsic differences in the polari- 
sation geometry between pulsars: thus, we have accepted at 
face value that the best determined PAo corresponds to the 
polarisation position angle at the minimum approach to the 
magnetic pole, which might not be the case, even for flat PA 
profiles, like that of PSR J2013+3845. 

4-2.2 The influence of subsets 

Arguably, the best, statistically significant observational 
support for pulsar spin-velocity alignment has come so far 
from the data of J05 — including the re-examination of those 
data by Rankin (2007). Follow up work, e.g. by Johnston et 
al. (2007), has provided a less convincing picture — although 
mitigating arguments were presented to explain this. Our 
statistical analysis included iVjos = 19 pulsars from the data 
set of J05, and the degree of observed alignment must there- 
fore have been influenced by the properties of that subset. 
In order to gauge how independent our analysis is compared 
to that of J05, we calculated the statistic D and it associ- 
ated probability p(D), as before, for data sets including a 
different number of pulsars from J05, ranging from to all 
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Figure 6. Probability of rejection of uniformity, under the KS test, as a function of (a) departure of PAo from the published values 
in units of the PA profile's peak-to-peak difference (as shown in Fig. EJl , and (b) fraction of values from Johnston et al. (2005) in the 
data. In both plots, the dashed black lines correspond to the ltr confidence interval around the probability of the mean KS statistic. The 
dotted grey lines mark the 68, 95.4 and 99.6% CLs of a normal distribution. 



19; in each of the data sets, we kept the same total number 
of pulsars, i.e. N — Njos = 35 pulsars, as a means to accom- 
modate the case where no J05 pulsars were included and 
to keep the samples statistically comparable. More specif- 
ically, the inclusion of different numbers of J05 pulsars in 
those data sets were performed in a bootstrapped manner, 
upon which a random sample (with replacement) is chosen 
to be part of the N pulsars in each data set. This process is 
re-iterated a large number of times by means of calculating 
the distribution of D for a particular ATjcm . As before, each 
bootstrapped data set was subjected to the same MC pro- 
cedure as for the original data set, in order to calculate the 
mean value of the KS statistic. The probability of rejecting 
uniformity, p(D), as a function of ATjob is plotted in Fig. [6)3. 
The error bars in that plot indicate the confidence in the cal- 
culated central value of p(D). Due to the positive skewness 
of the distribution of means from the MC, we chose to plot 
the probability of the median and its la asymmetric errors 
instead of that of the mean. 

As with the test on PAo, Fig. [6)3 shows a positive trend 
for increasing J05 fractions: i.e. from p(D) ~ 0.96, when no 
J05 value is included in the data, to p(D) > 0.99, when all 
J05 values are included. In addition, it is clear that the con- 
fidence interval of p(D) is significantly decreased towards 
higher values of Njos. Both these facts advocate the signif- 
icant impact of the J05 data set. It is true, however, that 
all bootstrapped data sets reject Ho at a higher that 2a 
level, even without including any J05 values — although, 
it should be noted that the confidence interval in the latter 
case extends to far below that level (i.e. p = 0.961q 13)- 



5 SUMMARY 

We have revisited the long-standing problem of pulsar spin- 
velocity alignment with at least twice as many pulsars as 
in other respective studies so far. A sample of 54 pulsars 
was selected from Parkes, Effelsberg and Lovell observations, 



with accurate PA profiles and reliable proper motions, which 
was used to examine the distribution of angles between the 
proper motion and spin-axis orientations. In addition, we 
included 3 pulsars for which the spin-axis orientation had 
been previously obtained by fitting the PWN torus sur- 
rounding the pulsar. Allowing for the possibility of orthogo- 
nal emission and the intrinsic 7t-ambiguity in measurements 
of the polarisation vector, we simulated the distribution of 
^ = PA„ — PAo mod 45° based on the measurement errors 
on the individual PAs. The distribution of it appears signif- 
icantly non-uniform, with 70% of the values residing below 
20° and the corresponding probability of rejecting unifor- 
mity under the KS test being 99.8%, which is a roughly 3<r 
rejection of the case that PA„ and PAo are uncorrelated. 
Ignoring the possibility of orthogonal emission on the other 
hand reveals that having spin and velocity orthogonal to 
each other is as likely as having them aligned. Notably, the 
probability distribution exhibits clear minima for offsets be- 
tween 30° and 60°, thus supporting the robustness of the 
alignment (or, indeed, orthogonality) between spin and ve- 
locity. 

The statistically well-defined shape of the distribution 
of ^P, shown as a histogram in Fig. [2] encouraged a compari- 
son between the observed distribution and a number of sim- 
ulated distributions generated from spin-aligned and spin- 
orthogonal velocity distributions, having a range of stan- 
dard deviations. The statistical comparison showed that the 
alignment scenario is more probable than spin-velocity or- 
thogonality, for the observed distribution of angles, although 
the confidence intervals in both cases are quite broad: a con- 
clusion that can already be drawn from the PDF of Fig. [2] 
Furthermore, this comparison allowed us to determine the 
spread that the kick-velocity distribution must have to be 
consistent with the observations: we derived a « 20° — 35° 
spread for the case of alignment. A much narrower, < 10° 
spread was found for the orthogonal case, which is consis- 
tent with previously published results that allow for binary 
kicks. 
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Finally, we checked the robustness of our conclusions 
against two systematic effects: we investigated (a) how toler- 
ant our conclusion is to departures from the best determined 
values of PA and (b) how much the inclusion of the John- 
ston et al. (2005) data set (the best statistically favourable 
evidence, yet) biases our conclusions. In both cases, we found 
that observed alignment is robust: excluding the J05 data set 
does little towards weakening the non-uniformity of the dis- 
tribution of although including it significantly improves 
the //o-rejection confidence; also, the observed alignment is 
favoured at > 95% confidence to within a 10% excursion 
from the published PAo values. It should be stressed that our 
robustness analysis could only probe systematics induced by 
human choice: systematics intrinsic to the pulsar emission, 
e.g. the polarisation orientation relative to the magnetic field 
lines, cannot be assessed without further information, such 
as X-ray imaging of PWN tori. 

The strong evidence for a correlation between spin and 
velocity vectors in our sample strengthens the conclusions of 
previous publications, but it also opens new questions about 
the physics behind the correlation. The signature of a kick 
process that favours alignment at pulsar birth should not 
be traceable for old pulsars (> 10, 000 kyr) whose veloci- 
ties have been stochastically altered by the Galactic grav- 
itational potential. On the other hand, a sample of young 
pulsars (< 1,000 kyr; like the ones selected in J05) should 
contain such a signature. Our sample did not make a distinc- 
tion between those age groups and yet a strong correlation 
was evident in the data (r c ~ 1 — 100, 000 kyr). The investi- 
gation of how robust the correlation is as a function of pulsar 
age is the subject of a follow-up study. 
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Figure 7. Polarisation and position-angle (PA) profiles for 51 pulsars taken from Johnston et al. (2005, 2007) and Carr (2007). The 
lower panel of each plot shows the total intensity profile (black) and the linear (red) and circular (blue) polarisation profiles. The top 
panels show the profile of the PA across the pulse, for all phase bins with S/N of the linear polarisation greater than 3, i.e. L/crj > 3. In 
each of the plots a vertical, dashed black line is drawn at PAo, i.e. the position angle of polarisation at the point of closest approach of 
the observer's LOS to the magnetic pole. Also shown in these plots is the position angle of the proper motion vector (dashed, green line). 
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Figure 2. Continued. 
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Figure 2. Continued. 
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